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Entanglement is an essential property of quantum many-body systems. However, its local detec¬ 
tion is challenging and was so far limited to spin degrees of freedom in ion chains. Here we measure 
entanglement between the spins of atoms located on two lattice sites in a one-dimensional Bose- 
Hubbard chain which features both local spin- and particle-number fluctuations. Starting with an 
initially localized spin impurity, we observe an outwards propagating entanglement wave and show 
quantitatively how entanglement in the spin sector rapidly decreases with increasing particle-number 
fluctuations in the chain. 

PACS numbers: 37.10.Jk, 67.85.-d, 75.10.Pq, 05.70.Ln 


Quantum many-body systems are distinct from their 
classical counterparts due to entanglement among their 
constituents [1, 2]. Especially in strongly correlated 
regimes, such as, in the vicinity of quantum phase tran¬ 
sitions [1, 3] or far away from equilibrium [4], the growth 
of entanglement with time or subsystem size seriously 
limits numerical simulations of complex quantum sys¬ 
tems [2]. Next to its importance on this fundamental 
level, entanglement is a valuable resource for quantum 
information and its microscopic control is required for 
most applications [5]. Experimentally, the measurement 
of entanglement is difficult given that full quantum state 
tomography requires extraordinary control and resources 
making it feasible only in small systems [6-8]. In larger 
or more complex many-body systems the, mere presence 
of entanglement can be inferred from macroscopic ob¬ 
servables, often relying on entanglement witnesses [1, 9]. 
Such a strategy has been applied for susceptibility mea¬ 
surements in solids [10], collective spin systems [11-16] 
or coupled superconducting qubits [17]. 

In Hubbard systems realized with ultracold atoms, en¬ 
tanglement in the on-site occupation number degree of 
freedom has been inferred from the visibility of a far- 
held interference pattern [18, 19]. A spatially resolved 
detection of entanglement has been recently proposed 
using the Renyi entropy [20, 21]. Extension of the Hub¬ 
bard model to two components introduces a spin degree 
of freedom [22, 23] such that spin-entanglement, in the 
sense of the concept of entanglement of particles [24, 25], 
can be present. First experiments with ultracold atoms 
showed that short-range coherent spin-dynamics can be 
controlled in bosonic [26-28] , as well as in fermionic sys¬ 
tems [29] . Spin exchange collisions in state selective opti¬ 
cal lattices have been used to realize collisional gates be¬ 
tween neighboring atoms [30] and global measurements 
indicated entanglement [31, 32]. However, a spatially 
resolved detection of either spin or occupation number 


entanglement in Hubbard models has still been an out¬ 
standing experimental challenge. 

Inspired by recent measurements in ion chains [8] , here 
we report on the spatially resolved detection of entan¬ 
glement among spin degrees of freedom in a two compo¬ 
nent Bose-Hubbard chain following a recent proposal [33] . 
Specifically, we used local detection [34, 35] and ma¬ 
nipulation [36] to study the dynamics of a single spin- 
impurity. Measurements in the longitudinal basis re¬ 
vealed the position of the impurity [37, 38], while trans¬ 
verse correlation measurements were used to extract its 
coherence. Combining both observables yields a lower 
bound for the concurrence [39, 40] in the spin-degrees of 
freedom of particles on two lattice sites [33] . We observe 
an outward propagating entanglement wave, the evolu¬ 
tion of which we follow up to a distance of six lattice sites. 
Importantly, the detected bound is valid despite on-site 
particle-number fluctuations and proves entanglement of 
particles [24, 25, 41] in our system [33]. To study the 
effect of occupation-number defects on the entanglement 
propagation in more detail, we developed a novel in-situ 
Stern-Gerlach imaging technique that yields information 
on both, the local spin and the occupation-number in one 
image. 

In our experiment, we realized ferromagnetic Heisen¬ 
berg spin chains with two-component ultracold bosonic 
atoms in an optical lattice [22, 23]. In the unity fill¬ 
ing Mott insulator regime and for equal inter- and intra¬ 
component scattering lengths, the system is described by 
the isotropic Heisenberg model H = — J ex 5f +1 + 

Si 1 + Si&i+i) + Hd, where the last term accounts 
for defects (holes or multiple occupancies) in the chain. 
The operators Sf denote the components of a spin-1/2 
operator at site i with a = x,y,z and J ex « 4 J 2 /U 
is the superexchange coupling for onsite interaction U 
and tunnel coupling J of the underlying Bose-Hubbard 
model. We prepared the system in a fully polarized state 
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FIG. 1. Single spin-impurity dynamics, (a) Schematics. The 
Bloch spheres represent the spin vector at different sites along 
the chain; their coupling with strength J ex is indicated by the 
arrows. The |t) (|D) states are shown as the gray areas at the 
north (south) poles. At time t — 0, a single spin-impurity is 
deterministically created in the center and subsequently prop¬ 
agates to both, left and right. For all times, the only non¬ 
vanishing first moment is (S'f), the longitudinal mean spin. 
The transverse direction is fully undefined (depicted for some 
later time t — r by the ring on the Bloch spheres), but the 
dynamics generates correlations that are represented by the 
changing color along the ring, where the same color indicates 
correlations in the spin direction between the sites. The size 
of the gray areas and the width of the ring are optimized for 
visibility and they have no physical meaning, (b) Representa¬ 
tive experimental measurements in the longitudinal (left) and 
transverse basis (right) after t = 35 ms. The Bloch sphere 
pictograms indicate the measurement direction. The images 
show exemplary single shot images. The histograms show the 
probability at each lattice site to find an atom in the respec¬ 
tive spin-resolved measurement. 


(i.e. the spin of all atoms points down) and created a sin¬ 
gle spin-up impurity on the central site. In this case, the 
interaction term SfSf +1 can be dropped and the Hamil¬ 
tonian reduces to the XX-model in the spin sector. The 
wavepacket of the spin-impurity dispersed with evolving 
time resulting in a build up and subsequent spreading 
of spin-entanglement along the chain [33, 42-44]. As an 
experimentally accessible entanglement measure, we use 
the concurrence between two lattice sites [39, 40], which 
measures entanglement on a scale between zero (no en¬ 
tanglement) and one (maximal entanglement). A con¬ 
venient lower bound, detectable with only global spin 
rotations, is given by = 2{2C it j - Pff) 1 / 2 ) [33]. 

The first term = (S^Sj-) = (( SfSJ) + (S- Sj))/2 
measures transverse (_L) spin correlations as the mean 
of the correlations in x and y directions. The latter 
term takes longitudinal spin correlations into account, 
where i s the joint probability to find spin-|D 


(ID) atoms at positions i and j. 

Similar to the experiments reported in [37, 38], 
we started with the preparation of a two-dimensional 
quantum-degenerate gas of typically 170 87 Rb atoms in 
the |D = |F s== 1, mp = — 1) state trapped in a single 
anti-node of a vertical (z-axis) optical lattice with a depth 
of V z = 20 F r , where E r = /i 2 /(8raa 2 at ) is the recoil 
energy with lattice spacing a\ at = 532 nm and atomic 
mass m. To prepare the system for the local spin flip, 
we adiabatically ramped up two horizontal lattices to 
Vx,y = 40 F r , which drove the gas deep into the Mott 
insulating phase [34]. At this stage, we extracted the 
temperature T = 0.08(3) U/k& from the density and dis¬ 
tribution of holes [34]. Notably, we improved our mini¬ 
mal temperature compared to earlier experiments by al¬ 
most a factor of two [37] resulting in a lower probability 
of 0.032(6) per site to be empty (in the central region 
of interest of nine sites). These empty sites might be 
either non- or doubly-occupied sites, which we cannot 
distinguish in our imaging procedure due to parity pro¬ 
jection [34]. A single line of atoms (in spatial ^/-direction) 
in the center of the Mott insulator was then transferred 
to the |D = \F = 2, rriF = —2) state using our single site 
addressing technique [36-38], thereby deterministically 
creating an impurity spin in each chain. The wavelength 
and polarization of the line-shaped addressing beam was 
adjusted such that the |D state was nearly unaffected 
while atoms in the |D state experienced an attractive 
potential. Next we decreased the lattice along the x- 
direction to V x = 10 E r within 50 ms while keeping the 
addressing beam on. The one-dimensional impurity dy¬ 
namics with J ex /Ti = 27r x 10 Hz, J /Ti = 27r x 39 Hz and 
U/h = 2tt x 800 Hz was then initiated by switching off 
the addressing beam in 1 ms. Here, J ex was directly ex¬ 
tracted from the impurity dynamics [37, 38], while J and 
U were calculated for the given lattice depths. Correc¬ 
tions due to density-induced tunneling resulted in an in¬ 
crease of J ex of about 20% compared to 4 J 2 /U [45]. 

After a variable evolution time £, we froze the atomic 
distribution by increasing the depth of all lattices. Se¬ 
lective imaging of the |D state was performed by in¬ 
verting the spin population using a global microwave 
sweep followed by a push out of the |2, —2) state on the 
cycling transition and subsequent site-resolved fluores¬ 
cence imaging of the remaining atoms [34]. The prob¬ 
ability P. 1 1| to find one atom at site i after this longi¬ 
tudinal (||) measurement corresponds to the probability 
for the atom to be in the |D state. Clear interference 
fringes with high visibility were observed in these mea¬ 
surements (see Fig. lb). In order to measure the anal¬ 
ogous quantity Pt_|_ in the transverse basis, we added 
a global 7t/2 rotation before the spin selective imaging. 
Spatially homogeneous magnetic field fluctuations ran¬ 
domize the transverse phase within less than 1 ms, re¬ 
sulting in rotational symmetry around the longitudinal 
spin axis. Consequently, the transverse spin distribution 
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FIG. 2. Transverse correlations, (a), (b) Experimental data 
of the correlations Cij and Cij for different evolution times. 
The strongest signal corresponds to the outward propagat¬ 
ing correlations between the sites d=z symmetrically located 
around the initial position (on the upper left to lower right 
diagonal), (c) Theoretical prediction for an ideal spin chain. 
Remarkable qualitative agreement between theory and exper¬ 
iment is visible in the spatial structure of the correlations, but 
the amplitude of the experimental signal is reduced (different 
color scales). 


is uniform across the chain with equal probabilities for 
the |t) and |t) states (Fig. lb). We ensured that the 
magnetic field homogeneity was better than 50mG/cm 
(0.2Hz/ai at ) [46] such that transverse spin correlations 
are preserved over experimental timescales of 100 ms. 
Each data point presented in this paper is extracted from 
typically 800 (1000) individual realizations of the spin 
chain in the longitudinal (transverse) case (with the ex¬ 
ception of the data shown in Fig. 4, where these numbers 
are five times lower). 

The spatially resolved measurement of transverse cor¬ 
relations Cij is the crucial step towards the detection 
of entanglement dynamics. Without defects in the spin 
chain, the operator S^~ directly relates to our experimen¬ 
tal observable = (S-~) + Hence, the transverse 
correlations Cij = P/j ± \ are given by the joint prob¬ 

ability Pj] a to find one atom at site i and one at site 
j in the transverse (a =_L) measurement. Imperfections 
will always decrease these detected correlations such that 
Cij provides a lower bound for them, even in an environ¬ 
ment of onsite atom number fluctuations [33]. Figure 2a 
shows the measured transverse correlations together with 
the theoretical prediction for the ideal XX-spin chain. A 
strong positive signal appears between sites +1 and — 1 
after an evolution time of 20 ms (1.2 6h/J ex ) and subse¬ 
quently these correlations spread further outwards. How¬ 
ever, compared to the ideal case a trend toward neg¬ 
ative values is visible even between far separated sites 
that should be uncorrelated given the short evolution 
times. A possible explanation for this lies in the non¬ 
perfect initial Mott insulators resulting in P/ ± < 0.5, 
which biases the measured Cij. This bias is removed 


in the modified transverse correlation Cij, defined as 
Cij = Pjj ± — Pi±P }_l [33]. In Fig. 2b, we show the 
measured Cij, now in remarkable agreement with the 
theoretical prediction, except for the smaller amplitude 
of the measured correlation signal. 

We now combine longitudinal and transverse correla¬ 
tion measurements to detect spin-entanglement in the 
system. This is achieved using a lower bound for the 
concurrence in the spin-1/2 degree of freedom and also 
for the entanglement of particles [24, 25, 41] obtained 
from [33] 

^j=2(2Cij~^pJpJ\). ( 1 ) 

Here, J*o is the joint probability of finding zero atoms 
on sites i and j in the longitudinal measurement. It has 
been shown in Ref. [33] that %j is a valid lower bound 
for the concurrence even in the case of fluctuating parti¬ 
cle numbers as long as the maximum on-site occupation 
number does not exceed two. This requirement is fulfilled 
in our experiments where the total atom number is tuned 
to yield a unity filled Mott insulator in the center of the 
trap. Assuming the worst case scenario, that the ob¬ 
served hole probability of 0.032(6) is only due to doubly 
occupied sites and an exponentially decreasing occupa¬ 
tion of higher excited states we expect a probability for 
triply occupied states of 10 -3 . With the weak additional 
assumptions of vanishing correlations between the site 
occupation numbers and between all degrees of freedom 
in the doubly occupied sector, a more efficient bound ^j 
on the concurrence can be obtained by replacing Cij by 
the modified transverse correlations Cij in Eq. (1) [33]. 
It is reasonable to assume that these conditions are ful¬ 
filled in the experiment for non nearest-neighbor sites and 
given the very low probability for double occupation. 

The results of the concurrence measurements are 
shown in Fig. 3a for pairs of sites symmetric around the 
initial impurity position. They reveal a buildup of en¬ 
tanglement in the spin chain leading to a peaking con¬ 
currence = 0.24(6) between the sites ±1 a\ at away 

from the center after 35 ms. For longer times, the con¬ 
currence peaks at larger distances showing an outward 
propagating entanglement wavefront. Using the bound 
^ij, we find finite entanglement up to distances of six 
lattice sites. Note that the bound for the concurrence is 
expected to be especially efficient for pairs of sites located 
symmetrically around the initial position, which is con¬ 
sistent with our observations (cf. Fig. 2). The observed 
concurrence closely resembles the transverse correlations 
that are shown in Fig. 3b for comparison. Its amplitude 
is only slightly decreased due to our finite fidelity in the 
preparation of the initial spin-impurity. Comparing the 
measured transverse correlations Cij quantitatively to 
the expectation for a perfect chain, we find good agree¬ 
ment after a constant down scaling of the theoretical cor¬ 
relation amplitude by 0.6. Such an effect has indeed been 
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FIG. 3. Propagation of an entanglement wave, (a) Experi¬ 
mental lower bounds for the concurrence (circles with 

solid line) and the more efficient bound (squares with 

dashed line) between the sites d =i versus time. From top to 
bottom i — 1 (red), 2 (blue), 3 (green) (b) Transverse spin 
correlations 4 C-i,+i for the same sites. Circles are the exper¬ 
imental data and the solid lines show the predictions for the 
defect-free ideal case, where 4(7-*,+* = ^-z,+z = The 

dashed lines are the ideal predictions scaled by 0.6. Error 
bars indicate the standard error of the mean (s.e.m.) 


predicted in an environment of holes and doubly occu¬ 
pied sites [33]. An decreasing entanglement signal with 
increasing site separation is also expected in the ideal 
case due to the dispersion of the single spin-impurity. 

The natural question arising from these results is the 
impact of defects in the chain on coherence and entan¬ 
glement in the spin sector. To address this question, 
we refined our measurement technique to simultaneously 
detect the local occupation number in both spin states, 
which also directly gives access to hole defects in the sys¬ 
tem. Our technique is based on an in-situ Stern-Gerlach 
like measurement [47], which spatially separates the |t) 
and \l) spins prior to the detection [33]. Since imaging 
of the atoms is restricted to a single plane, this required 
the preparation of single isolated one-dimensional (Id) 
systems such that the Stern-Gerlach separation can be 
done transversally. The experimental sequence closely 
followed the procedure described earlier, but we addition¬ 
ally use the local addressing system to remove all atoms 
but those in a single Id tube prior to the preparation of 
the spin-impurity. After the spin evolution, we switched 
off the lattice in ^/-direction (transverse to the chain) and 
applied a magnetic field gradient along the same direc¬ 
tion. Due to the different magnetic moments of the |t) 
and 1 1) states, the two components were spatially sepa¬ 
rated. Next, we switched the lattice in ^/-direction back 
on in 75 ms to 10 E r . Finally, we used our standard spin- 
insensitive fluorescence detection to observe the atoms 
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FIG. 4. Impact of defects on spin-entanglement. Compari¬ 
son of images of single spin chains taken with in-situ Stern- 
Gerlach imaging of both spin states (a) to the standard spin- 
resolved imaging in which one spin state was removed before 
detection (b). In the former case the magnetic field gradient 
was in vertical direction, pushing the 17}-spins down and the 
1 1 )-spins up from the initial position (white line), (c) Lower 
concurrence bound ^_i 5 +i after 35 ms evolution for the full 
dataset (point in the left gray part) and for subsamples postse- 
lected to 0, 1, 2 and 3 holes in the chain. The red dashed line 
indicates the value expected for the perfect XX-spin chain, 
the black dashed line is the boundary for the presence of en¬ 
tanglement. The inset schematically shows the Stern-Gerlach 
separation of the two spin states in the magnetic field gradient 
B' (depicted by the left triangle). Error bars 1 s.e.m. 


locally. In the obtained image, the position along x de¬ 
fines the position in the spin chain and the position in 
^/-direction identifies the spin state (Fig. 4a). This mea¬ 
surement is very challenging, both, due to the increased 
complexity of the preparation and due to the lower statis¬ 
tics; we measure only a single chain per experimental run 
as compared to nine chains in the simpler standard pro¬ 
tocol. Therefore, we limited ourselves to a fixed evolu¬ 
tion time of 35 ms, the setting for which the concurrence 
peaked between the ±1 sites. 

The results for the concurrence extracted from the 
spin-resolved measurements are shown in Fig. 4c. For 
the analysis, we discarded all pictures with more than one 
atom per 7/-tube (approx. 25%), as those were predomi¬ 
nantly caused by imperfect preparation of the single Id 
system. After this, the concurrence signal agrees 

with the previously measured one within experimental 
uncertainty - an indication that double occupancies had 
little impact on the measurements reported above. Next, 
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we further postselected the data to a fixed number of 
holes (between zero and three) within the central nine 
sites. For the zero hole data set, our measured concur¬ 
rence indeed matches the expectation for an ideal XX- 
spin chain while it rapidly decreases for increasing hole 
number. For three holes, i.e. a hole density of 30%, no 
entanglement is detected any more showing that defects 
critically affect the coherence in the spin sector. 

In conclusion, we have experimentally measured bipar¬ 
tite entanglement in the spin degree of freedom between 
two sites of an optical lattice filled with ultracold atoms. 
Using a novel detection method we have shown that atom 
number defects are critical for the transverse coherence 
and hence for entanglement in the spin sector. This 
method also gives access to inter-spin correlations, an im¬ 
portant quantity to characterize two-component bosons 
in the superfluid regime [37]. Our measurements pave the 
way toward in-depth studies of entanglement in quantum 
many-body entanglement systems [1] and mark the first 
steps toward controlled entanglement transfer across spin 
wires [42, 43, 48]. 
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MPG, EU (UQUAM) and the Korber Foundation. 
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